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Effects  of  Concave  Curvature  on 
Boundary  Layer  Transition  Under 
High  Freestream  Turbulence 
Conditions 

An  experimental  investigation  has  been  earned  out  on  a  transitional  boundary  layer 
subject  to  high  (initially  9%)  freestream  turbulence,  strong  acceleration  (K  =  (v/U(J 
X(dUw/dx)  as  high  as  9 X10~6),  and  strong  concave  curvature  (boundary  layer  thick¬ 
ness  between  2%  and  5%  of  the  wall  radius  of  curvature).  Mean  and  fluctuating  velocity 
as  well  as  turbulent  shear  stress  are  documented  and  compared  to  results  from  equivalent 
cases  on  a  flat  wall  and  a  wall  with  milder  concave  curvature.  The  data  show  that 
curvature  does  have  a  significant  effect,  moving  the  transition  location  upstream,  increas¬ 
ing  turbulent  transport,  and  causing  skin  friction  to  rise  by  as  much  as  40%.  Conditional 
sampling  results  are  presented  which  show  that  the  curvature  effect  is  present  in  both  the 
turbulent  and  nonturbulent  zones  of  the  transitional  flow.  [DOI:  10.1115/1.1522410] 


Introduction 

Boundary  layer  transition  is  influenced  by  a  number  of  factors, 
including  streamwise  pressure  gradient,  freestream  turbulence  in¬ 
tensity  (FSTI),  surface  curvature,  compressibility,  and  surface 
roughness.  In  gas  turbine  environments,  Mayle  [1]  states  that  a 
substantial  fraction  of  the  boundary  layer  on  both  sides  of  a  gas 
turbine  airfoil  may  be  transitional.  The  extended  transition  zones 
exist  due  to  strong  favorable  pressure  gradients,  found  on  both  the 
pressure  side  and  the  leading  section  of  the  suction  side,  which 
stabilize  the  boundary  layer  and  delay  transition  in  spite  of  the 
high  freestream  disturbance  levels  in  gas  turbine  environments. 
The  ability  to  model  and  predict  high  FSTI  transition  is  important 
since  heat  transfer  rates,  skin  friction  coefficients,  and  in  some 
cases  boundary  layer  separation  depend  strongly  on  the  state  of 
the  boundary  layer  with  respect  to  transition.  Improved  transition 
models  and  turbine  designs  depend,  therefore,  on  a  better  under¬ 
standing  of  high  FSTI  transition. 

The  role  of  surface  curvature  on  high  FSTI  transition  is  not 
fully  understood.  Mayle  [1]  states  that  transition  onset  is  con¬ 
trolled  by  the  freestream  turbulence  and  the  periodic  unsteadiness 
caused  by  wakes  from  upstream  airfoils.  He  notes  that  the  turbu¬ 
lent  spot  production  rate,  which  determines  the  length  of  the  tran¬ 
sition  region,  is  controlled  primarily  by  the  freestream  disturbance 
and  the  streamwise  pressure  gradient,  but  that  curvature  may  play 
a  secondary  role.  The  extent  of  this  secondary  role  is  largely  un¬ 
known,  as  only  limited  data  exist  from  convex  surfaces  and  “vir¬ 
tually  no  reliable  data”  are  available  from  concave  surfaces  for 
computation  of  spot  production  rates,  [1],  Determining  the  signifi¬ 
cance  of  curvature  is  potentially  important  for  improved  turbine 
design.  If  curvature  is  not  significant,  both  modeling  of  transi¬ 
tional  boundary  layers  and  future  experimental  studies  may  be 
simplified,  as  results  from  flat-plate  studies  may  be  directly  ap¬ 
plied  to  curved  airfoils.  If  curvature  is  found  to  be  significant,  its 
effects  should  be  incorporated  into  future  transition  models. 

Reasonable  arguments  can  be  made  both  for  and  against  the 
importance  of  curvature  on  high  FSTI  transition.  Low  FSTI  re¬ 
sults  suggest  that  curvature  is  important.  Gortler  [2]  determined 
theoretically  that  convex  curvature  is  stabilizing  and  that  concave 
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curvature  is  destabilizing.  He  predicted  the  formation  of  stream- 
wise  vortices  (now  known  as  Gortler  vortices)  on  concave  sur¬ 
faces.  Liepmann  [3]  confirmed  these  results  experimentally  under 
low  FSTI  conditions.  He  showed  that  convex  curvature  only 
slightly  delays  transition,  but  that  concave  curvature  can  cause 
transition  to  occur  significantly  earlier.  Transition  occurred  when 
the  Gortler  number,  G,  was  between  6  and  9,  decreasing  with  the 
FSTI.  These  results  have  been  confirmed  in  several  subsequent 
studies.  Floryan  [4]  and  Saric  [5]  provide  reviews.  Volino  and 
Simon  [6]  provide  recent  documentation  from  a  zero  pressure  gra¬ 
dient  case,  and  Finnis  and  Brown  [7]  considered  a  favorable  pres¬ 
sure  gradient  case.  In  a  low  FSTI,  fully  turbulent  boundary  layer, 
Simonich  and  Moffat  [8]  showed  that  concave  curvature  of 
strength  comparable  to  gas  turbine  conditions  resulted  in  about  a 
20%  enhancement  in  heat  transfer  over  flat-plate  results. 

While  the  low  FSTI  studies  show  the  importance  of  concave 
curvature,  the  results  of  high  FSTI  investigations  are  less  clear. 
Data  from  fully  turbulent  boundary  layers  again  suggest  that  cur¬ 
vature  is  important.  Kim  et  al.  [9]  documented  cases  with  inlet 
FSTI  of  8%  on  both  flat  and  concave  walls,  and  showed  increased 
turbulent  activity  and  an  enhancement  of  heat  transfer  by  about 
16%  on  the  concave  wall.  Kestoras  and  Simon  [10]  documented  a 
boundary  layer  with  8%  inlet  FSTI  moving  from  a  concave  wall 
onto  a  flat  wall.  They  showed  an  almost  immediate  drop  in  turbu¬ 
lence  within  the  boundary  layer  as  the  flow  moved  onto  the  flat 
wall.  They  postulated  a  synergistic  effect  between  curvature  and 
freestream  turbulence,  stating  that  when  concave  curvature  and 
high  FSTI  are  combined,  the  freestream  eddies  are  able  to  pen¬ 
etrate  closer  to  the  wall  than  in  a  flat  wall  flow,  resulting  in  sig¬ 
nificantly  higher  transport  within  the  boundary  layer. 

Only  a  few  studies  have  documented  transition  with  high  FSTI. 
On  surfaces  subject  to  zero  streamwise  pressure  gradients,  Blair 
[11],  Sohn  and  Reshotko  [12],  and  Kim  et  al.  [9]  all  showed  that 
at  FSTI  above  about  3%,  transition  occurred  rapidly  near  the  lead¬ 
ing  edge  of  a  test  surface.  Hence,  little  documentation  of  transition 
is  provided  in  these  cases.  Kim  et  al.  [9]  included  a  concave  cur¬ 
vature  case,  and  found  no  evidence  of  Gortler  vortices.  They  pro¬ 
posed  that  the  enhanced  turbulent  transport  caused  by  the  high 
FSTI  and  early  transition  lowered  the  effective  Gortler  number 
and  suppressed  the  formation  of  the  vortices.  Since  Gortler  vorti¬ 
ces  change  the  transition  mechanism  and  promote  early  transition 
under  low  FSTI  conditions,  [6],  their  absence  under  high  FSTI 
might  suggest  that  curvature  does  not  play  a  strong  role  in  high 
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FSTI  transition.  Riley  et  al.  [13]  considered  zero  pressure  gradient 
cases  on  concave  walls  with  inlet  FSTI  as  high  as  7.2%.  Their 
results  indicate  that  with  very  strong  curvature  ((r/ ds)°'5< 20), 
transition  start  is  delayed  relative  to  flat  wall  results.  They  re¬ 
ported  the  presence  of  Gortler  vortices,  and  attributed  the  delayed 
transition  to  increased  near-wall  velocity  gradients  caused  by  the 
vortices.  In  cases  with  curvature  representative  of  gas  turbine  con¬ 
ditions  ((r/0s)a5> 30),  as  indicated  by  Mayle  [1]),  Riley  et  al. 
[13]  showed  no  significant  difference  in  transition  start  between 
concave  and  flat  wall  results.  They  did  not  provide  a  correlation 
for  transition  end. 

Extended  transition  regions  do  occur  when  high  FSTI  is  com¬ 
bined  with  strong  acceleration.  Results  are  available  from  turbine 
cascade  and  rotating  rig  studies  (e.g.,  Halstead  et  al.  [14]),  but  in 
these  cases  the  effects  of  airfoil  curvature,  while  present,  cannot 
be  isolated.  Volino  and  Simon  [15-17]  considered  transition  along 
a  concave  wall  with  inlet  FSTI  of  8%  and  acceleration  with  K  as 
high  as  9X  1CT6.  Acceleration  rates,  Reynolds  numbers  and  FSTI 
were  typical  of  the  pressure  side  of  a  gas  turbine  airfoil.  The 
strength  of  curvature,  (r/0s)°'5  =  53,  was  mild  relative  to  gas  tur¬ 
bine  conditions.  An  extended  transition  region,  with  intermittent 
turbulent  and  nonturbulent  zones,  covered  most  of  the  test  surface. 
Although  detailed  documentation  of  the  transition  was  provided, 
no  comparison  case  from  a  flat  wall  was  available,  so  the  effect  of 
the  concave  curvature  could  not  be  determined.  Volino  and  Simon 
[17]  stated  that  the  curvature  effect  was  probably  small,  since  the 
strong  acceleration,  by  suppressing  the  growth  of  the  boundary 
layer,  kept  the  strength  of  curvature  low.  In  contrast,  Volino  [18] 
presented  a  model  incorporating  freestream  turbulence  and  curva¬ 
ture  effects  and  used  it  to  predict  that  the  curvature  effect  in 
the  Volino  and  Simon  [15-17]  case  might  be  substantial.  As 
stated  above,  no  experimental  data  were  available  to  verify  the 
prediction. 

The  present  study  addresses  the  significance  of  concave  curva¬ 
ture  directly.  The  case  presented  by  Volino  and  Simon  [15-17] 
has  been  reproduced  on  a  flat  test  wall,  and  a  concave  wall  with 
strong  curvature  ((r/ds)05=  27).  The  flat  wall  results  were  pre¬ 
sented  in  detail  by  Volino  et  al.  [19].  In  the  present  paper  the 
concave  wall  results  are  presented  and  compared  to  the  flat  wall 
results  and  the  weaker  concave  curvature  case  of  Volino  and  Si¬ 
mon  [15-17], 

Experiments 

Facility  and  Measurements.  All  experiments  were  con¬ 
ducted  in  the  same  low  speed  wind  tunnel  described  in  detail  by 
Volino  et  al.  [19]  and  shown  in  Fig.  1.  The  only  change  to  the 
facility  for  the  present  study  was  the  replacement  of  a  flat-plate 
test  section  with  a  curved,  converging  channel.  The  facility  is  very 
similar  to  that  used  by  Volino  and  Simon  [17],  again  with  the 
exception  of  the  curvature  of  the  test  wall.  One  side  of  the  present 
test  section  is  a  concave  curved  Plexiglas  surface  of  0.69  m  width 
and  0.80  m  length,  which  serves  as  the  test  wall.  This  wall  has  a 
constant  radius  of  curvature  of  25.4  cm.  Pressure  taps  are  installed 
along  its  span  wise  centerline.  At  the  leading  edge  of  the  test  wall 
a  slot  is  used  to  bleed  off  the  boundary  layer  which  grows  in  the 
development  section  between  the  turbulence  generator  and  the 
leading  edge.  Opposite  the  test  wall  is  a  flexible  convex  wall  that 
can  be  adjusted  to  set  the  desired  pressure  gradient  along  the  test 


Fig.  1  Schematic  of  the  test  facility 
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wall.  For  the  present  study  the  inlet  velocity  is  set  to  4.6  m/s  and 
the  velocity  gradient  along  the  wall  is  held  constant  at  13.9  s  , 
matching  the  previous  flat  wall,  [19],  and  mild  concave  curvature, 
[17],  cases.1  The  acceleration  parameter,  K,  drops  from  a  maxi¬ 
mum  of  9X  1CT6  at  the  inlet  to  the  test  section  to  1  X  10~6  at  the 
last  measurement  station.  The  FSTI  at  the  inlet  to  the  test  section 
is  8.6%,  based  on  all  three  components  of  the  fluctuating  velocity. 
Further  details  of  the  freestream  turbulence  including  spectra  and 
length  scales  are  given  in  Volino  et  al.  [19].  Within  the  test  sec¬ 
tion,  the  streamwise  component  of  the  freestream  turbulence,  u'x  , 
remains  nearly  constant  at  ~0.35  m/s,  but  vL  increases  from  0.50 
m/s  to  0.73  m/s.  On  the  flat  wall,  decreased  in  the  streamwise 
direction  while  v ^  remained  nearly  constant.  It  should  be  noted, 
however,  that  these  values  are  based  on  data  at  the  measurement 
point  most  distant  from  the  wall  (y  =  30  mm) .  While  n'„  and  v  4 
on  the  flat  wall  and  ii'x  on  the  concave  wall  reach  their  freestream 
value  by  this  location,  v ^  is  still  changing  with  distance  from  the 
wall  at  several  streamwise  stations  on  the  concave  wall. 
Freestream  turbulence  intensity  (normalized  using  the  local 
freestream  velocity)  drops  in  the  streamwise  direction  in  both  the 
concave  and  flat  wall  cases,  mainly  due  to  the  increasing 
freestream  velocity.  Values  of  K  through  the  test  section,  measure¬ 
ment  locations,  and  other  parameters  are  shown  in  Table  1. 

Velocity  profile  measurements  were  made  at  nine  streamwise 
stations  along  the  centerline  of  the  test  wall  using  a  single  sensor, 
boundary  layer  type  hot-wire  probe  (TSI  model  1218-T1.5),  a 
boundary  layer  cross-wire  probe  (TSI  model  1243-T1.5)  and  a 
constant-temperature  hot-wire  anemometer  (TSI  model  IFA-100). 
The  probes  were  moved  normal  to  the  wall  at  each  station  using  a 
manual  traversing  stage  with  a  minimum  increment  of  10  /rm.  At 
each  position  in  the  velocity  profiles,  data  were  acquired  for  26 
seconds  at  a  20  kHz  sampling  rate  (2 19  data  points).  The  hot-wire 
signals  were  low-pass  filtered  at  10  kHz.  All  raw  data  were  stored. 
Sampling  at  20  kHz  provided  essentially  continuous  velocity 
traces  for  subsequent  processing  as  detailed  in  Volino  et  al.  [19], 
Mean  and  rms  fluctuating  velocities  in  the  streamwise  and  wall- 
normal  directions  were  computed  from  the  instantaneous  data. 
Uncertainties  in  these  quantities  are  3-5%  except  in  the  very  near 
wall  region  (y  +  <5)  where  near-wall  corrections  (Wills  [22])  are 
applied  to  the  mean  velocity.  Uncertainty  in  the  turbulent  shear 
stress,  —u'v',  is  10%.  Skin  friction  coefficients  were  determined 
using  a  technique  whereby  t„  was  adjusted  until  the  U+  versus 
y+  data  fit  profiles  computed  using  near-wall  similarity,  as  de¬ 
scribed  by  Volino  and  Simon  [23],  This  technique  accounts  for 
pressure  gradient  effects  on  the  near  wall  profile.  Uncertainty  in 
Cf  is  8%.  Boundary  layer  thicknesses  were  determined  from  the 
mean  velocity  profiles.  Uncertainties  in  the  momentum  and  dis¬ 
placement  thicknesses  are  10%.  These  uncertainties  include  bias 
errors  which  tend  to  cancel  such  that  the  uncertainty  in  the  shape 
factor,  H,  is  7%.  The  intermittency  function,  which  indicates 
whether  the  boundary  layer  is  instantaneously  turbulent  or  nontur¬ 
bulent,  was  calculated  using  both  the  instantaneous  streamwise 
velocity,  u,  and  the  instantaneous  turbulent  shear  stress,  —  u'v' , 
as  described  by  Volino  et  al.  [19].  The  overall  uncertainty  in  the 
time-averaged  intermittency,  y,  is  0. 1 .  The  instantaneous  intermit¬ 
tency  function  was  used  to  separate  the  turbulent  and  nonturbulent 


'The  nominal  “freestream”  velocity  for  the  concave  wall  cases  is  not  so  clear  as 
in  the  corresponding  flat  wall  case.  On  a  flat  wall,  the  velocity  gradient  dUldy 
approaches  zero  outside  the  boundary  layer,  so  the  freestream  velocity  is  independent 
of  the  distance  from  the  wall.  Curvature  induces  a  nonzero  dU/dy,  which  under  low 
FSTI  conditions  approaches  the  potential  flow  solution  for  a  curved  channel.  With 
high  FSTI,  however,  the  combination  of  streamwise  curvature  and  freestream  turbu¬ 
lence  results  in  cross  transport  of  momentum,  as  described  by  Eckert  [20],  and 
deviation  from  the  potential  flow  solution.  The  gradient  dU/dy  is  approximately 
constant  in  the  freestream,  so  following  the  procedure  of  Kestoras  and  Simon  [21],  a 
straight  line  is  fit  to  the  mean  velocity  data  in  the  freestream  at  each  station  and 
extrapolated  to  the  wall  (y  =  0).  This  extrapolated  velocity,  designated  Uw  ,  is  used 
as  the  freestream  velocity  for  the  concave  wall  case. 
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Table  1  Flow  parameters  at  measurement  stations 
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Fig.  2  Intermittency  profiles  based  on  u 


Fig.  3  Intermittency  factor  versus  streamwise  distance 


zones  of  the  boundary  layer  for  conditional  sampling.  In  the  fig¬ 
ures  which  follow,  data  points  are  shown  for  the  nonturbulent 
zone  only  when  the  local  y<95%,  and  for  the  turbulent  zone  only 
when  y>5%.  Bulk  parameters  such  as  C,  are  presented  when 
ypk< 95%  and  ypk>5%  for  the  nonturbulent  and  turbulent  zones, 
respectively. 


Results 

Intermittency.  Intermittency  profiles  for  the  nine  measure¬ 
ment  stations  are  shown  in  Fig.  2.  These  results  were  computed 
using  the  instantaneous  streamwise  velocity,  u.  Although  not  pre¬ 
sented  here,  intermittency  profiles  based  on  the  instantaneous  tur- 


Table  2  Transition  start  and  end  locations 


Case 

( r/ef 5 

x,  [m] 

xe  [m] 

ncrxlO11 

Flat  [19] 

00 

0.29 

0.98 

4.2 

Mild  Curve  [17] 

53 

0.23 

0.90 

5.5 

Strong  Curve 

27 

0.12 

0.75 

10. 
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bulent  shear  stress  always  agreed  with  those  in  Fig.  2  to  within 
0.10  and  were  generally  within  0.05.  The  intermittency  remains 
nonzero  well  away  from  the  wall,  to  yl  S*  as  high  as  40.  In  the  flat 
wall  case,  Volino  et  al.  [19]  showed  y  approaching  zero  by 
y/<5*=10.  Something  in  the  concave  wall  case  promotes  turbu¬ 
lence  well  away  from  the  wall.  The  peak  intermittency  in  each 
profile  is  shown  in  Fig.  3  and  compared  to  results  from  the  flat 
wall  and  mild  curvature  cases.  The  intermittency  remains  fairly 
low  (ypk~10%)  for  the  first  two  stations  of  the  strong  curvature 
case  and  then  starts  to  rise.  Note  that  even  at  the  first  station,  the 
intermittency  is  higher  than  that  measured  in  the  other  two  cases. 
For  all  these  cases,  the  beginning  of  the  rise  in  y  corresponds  to  K 
dropping  below  3X1CT6.  In  low  FSTI  boundary  layers,  K> 3 
X  1CF6  leads  to  relaminarization  (lones  and  Launder  [24]).  By  the 
last  station,  transition  is  complete  in  the  strong  curvature  case, 
whereas  the  flat  and  mildly  curved  wall  cases  are  only  approach¬ 
ing  the  end  of  transition.  Following  the  technique  of  Narasimha 
[25],  as  modified  by  Volino  and  Simon  [26],  the  function 

f(ypk)=(-Hi-yPk))m  (1) 

can  be  computed  based  on  the  peak  intermittency  at  each  station 
and  plotted  versus  streamwise  location.  The  data  in  these  coordi¬ 
nates  tend  to  lie  along  a  straight  line.  The  line  may  be  extrapolated 
t°  /(yp<t=0)  =  0  and /(ypt  =  0.99)  =  2.146,  corresponding  to  the 
beginning  and  end  of  transition.  Transition  start  and  end  locations 
for  the  three  cases  are  listed  in  Table  2.  Figure  4  shows  ypk  plotted 
versus  dimensionless  streamwise  location  within  transition.  Also 
shown  is  the  theoretical  curve 

rM  =  exp(— 4.6(^)  )  (2) 

based  on  the  Dhawan  and  Narasimha  [27]  transition  model. 
Agreement  between  the  data  and  theory  is  good.  Differences  from 
Eq.  (2)  are  seen  at  low  y  in  all  the  cases.  This  is  termed  “subtran¬ 
sition”  [25]  and  is  expected  in  favorable  pressure  gradient  flows. 
The  Dhawan  and  Narasimha  [27]  model  assumes  a  concentrated 
breakdown  or  formation  of  turbulent  spots  at  the  transition  start 
location.  In  fact,  spots  are  formed  over  some  finite  distance  in  the 
streamwise  direction  in  a  distributed  breakdown.  Favorable  pres¬ 
sure  gradients  stabilize  the  boundary  layer  and  extend  this  region, 
resulting  in  the  subtransition.  The  presence  of  the  subtransition 
means  that  the  actual  start  of  transition  does  not  occur  at  the  xs 
location  given  in  Table  2.  The  xs  and  xe  value  do,  however,  indi¬ 
cate  the  rate  at  which  transition  proceeds.  The  mild  curvature  case 
completes  transition  slightly  upstream  of  the  flat  wall  case,  and 
the  strong  curvature  case  completes  transition  well  upstream  of 


Fig.  4  Peak  intermittency  in  profile  versus  dimensionless 
streamwise  location 

Journal  of  Fluids  Engineering 


Fig.  5  Mean  velocity  profiles  in  wall  coordinates:  (a)  compos¬ 
ite;  (b)  nonturbulent;  (c)  turbulent 


the  other  two  cases.  Since  curvature  is  the  only  difference  between 
the  cases,  it  is  clear  that  concave  curvature  has  a  significant  effect 
on  transition. 

The  dimensionless  turbulent  spot  propagation  rate  (Mayle  [1]) 
may  be  computed  as 

A.6v2U„ 
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Fig.  6  Mean  velocity  profile  for  Station  5 


where  is  the  average  freestream  velocity  in  the  transition  re¬ 
gion.  In  the  strong  curvature  case  fit r=1.0Xl0-1°.  For  the  flat 
wall  case  and  the  mild  curvature  case,  na  was  4.2X  1CT11  and 
5.5 X  10_  ,  respectively.  These  values  are  in  agreement  with  the 
trends  reported  by  Mayle  [1]  for  other  favorable  pressure  gradient 
cases  and  are  three  to  five  times  less  than  would  be  expected  for  a 
zero  pressure  gradient  case  with  the  same  FSTI  at  the  beginning 
of  transition.  It  appears  that  wall  curvature,  FSTI  and  pressure 
gradient  are  all  important  parameters  in  determining  the  spot 
propagation  rate. 

Mean  Velocity  Profiles.  Mean  velocity  profiles  for  the  nine 
measurement  stations  of  the  strong  curvature  case  are  presented  in 
Fig.  5  in  wall  coordinates.  Figure  5(a)  shows  the  composite  (un¬ 
conditioned)  profiles.  Figures  5  (b)  and  5(c)  show  turbulent  and 
nonturbulent  zone  results  from  conditional  sampling  of  the  data. 
The  composite  profiles  vary  only  slightly  through  Station  5.  For 
the  remaining  downstream  stations,  the  profiles  collapse  to  a 
turbulent-like  shape.  Even  the  most  upstream  stations  do  not  ex¬ 
hibit  the  laminar-like  shape  that  was  seen  in  the  flat  wall  results  of 
Volino  et  al.  [19].  At  Stations  7-9,  where  the  acceleration  rate  has 
dropped  to  K<  1.5X  10~6  and  transition  is  complete  or  nearly  so, 
the  profiles  agree  with  the  zero  pressure  gradient  law  of  the  wall, 
although  the  log-law  region  is  quite  small.  The  short  log-law  re¬ 
gion  is  a  consequence  of  the  combined  concave  curvature  and 
high  FSTI  that  gives  rise  to  significant  momentum  transport  and 
result  in  a  large,  negative  wake  strength.  The  nonturbulent  zone 
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Fig.  8  Shape  factor  versus  streamwise  distance 


profiles  (Fig.  5(b))  exhibit  a  laminar-like  shape,  even  at  the  end 
of  transition,  as  was  also  seen  in  the  flat  wall  case.  The  turbulent 
zone  profiles  (Fig.  5(c))  exhibit  a  more  turbulent-like  shape,  and 
the  last  five  stations  collapse.  Note  that  the  data  at  even  the  most 
downstream  station  falls  slightly  below  the  zero  pressure  gradient 
law  of  the  wall,  again  due  to  the  enhanced  momentum  transport  in 
the  outer  part  of  the  boundary  layer.  Examination  of  the  results  in 
[19]  shows  that  this  effect  is  not  as  strong  on  the  flat  wall,  where 
there  is  better  agreement  with  the  zero-pressure  gradient  law  of 
the  wall. 

Figure  6  shows  the  differences  between  the  composite,  nontur¬ 
bulent,  and  turbulent  profiles  at  Station  5,  in  the  center  of  the 
transition  zone  (ypk  =  0.67).  Also  shown  for  comparison  are  the 
results  for  the  flat  wall  at  the  same  streamwise  position,  which 
was  closer  to  the  start  of  transition  (ypk  —  0.\l).  On  the  concave 
wall,  velocities  in  the  near  wall  region  are  slightly  higher  in  the 
turbulent  zone  than  in  the  nonturbulent  zone  due  to  higher  levels 
of  turbulent  mixing.  Differences  between  the  zones  are  not  as 
dramatic  as  on  the  flat  wall,  however.  While  this  might  be  attrib¬ 
uted  to  the  difference  in  the  intermittency  for  the  two  cases  at  this 
streamwise  position,  it  is  also  observed  when  comparing  profiles 
from  the  two  cases  with  the  same  intermittency  and  different 
streamwise  locations.  The  enhanced  mixing  induced  by  the  curva¬ 
ture  results  in  higher  near  wall  velocities,  particularly  in  the  non¬ 
turbulent  zone.  In  the  concave  wall  case,  as  in  the  flat  wall  case, 
the  nonturbulent  zone  of  the  boundary  layer  never  behaves  as  if  it 
were  laminar  due  to  the  freestream  buffeting  the  boundary  layer. 
This  is  explained  further  in  Volino  [18]. 


Re, 


Fig.  9  Comparison  of  skin-friction  coefficient 
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Fig.  10  Fluctuating  streamwise  velocity  profiles  in  wall  coor¬ 
dinates:  (a)  composite;  (ft)  nonturbulent;  (c)  turbulent 


Boundary  Layer  Growth.  Figure  7  shows  the  momentum 
thickness  plotted  versus  streamwise  position.  On  the  concave 
wall,  the  composite  flow  momentum  thickness  remains  nearly 
constant  through  the  first  five  stations  and  then  increases  as  K 
drops  and  transition  nears  completion.  In  the  nonturbulent  zone, 
the  momentum  thickness  drops  somewhat  in  the  streamwise  direc¬ 
tion.  In  contrast,  the  turbulent  zone  momentum  thickness  in¬ 
creases  continuously.  This  is  likely  due  to  turbulent  entrainment  at 
the  boundary  layer  edge,  which  appears  to  be  more  significant  on 
the  concave  wall  than  the  flat  wall. 
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Fig.  11  Fluctuating  streamwise  velocity  profile  for  Station  5 


The  shape  factor,  H,  shown  in  Fig.  8,  is  an  indicator  of  the  state 
of  the  boundary  layer  with  respect  to  transition.  The  composite 
flow  shape  factor  drops  from  about  1.8  at  Station  1  to  1.45  at 
Station  9.  The  decrease  in  H  tracks  the  increase  in  ypk  (Fig.  3).  In 
the  nonturbulent  zone  H  increases  slightly  from  1.8  to  about  2.0. 
A  low  FSTI  laminar  boundary  layer  on  a  flat  wall  with  the  same 
pressure  gradient  would  have  a  shape  factor  of  about  2.4.  The 
shape  factors  below  the  laminar  value  agree  with  Fig.  6,  which 
shows  that  the  high  FSTI  makes  the  nonturbulent  zone  velocity 
profiles  appear  less  laminar-like  and  more  like  those  in  the  turbu¬ 
lent  zone.  In  the  turbulent  zone,  H  drops  from  about  1.5  to  1.3. 
These  values  are  significantly  less  than  those  for  the  flat  wall  case 
and  suggest  that  the  wall  curvature  promotes  greater  momentum 
transport  in  the  turbulent  zone  even  in  a  high  FSTI,  strongly  ac¬ 
celerated  boundary  layer. 

Skin  Friction  Coefficients.  Skin  friction  coefficients,  Cf, 
were  computed  from  the  mean  velocity  profiles  and  are  plotted 
versus  Re,  in  Fig.  9.  Composite  flow  results  are  presented  for  the 
flat  wall,  mild  and  strong  curvature  cases,  and  conditional  sam¬ 
pling  results  are  presented  for  the  flat  wall  and  strong  curvature 
cases.  The  striking  characteristic  of  the  figure  is  that  the  nontur¬ 
bulent,  turbulent,  and  composite  skin  friction  values  in  the  strong 
curvature  case  are  all  significantly  higher  than  in  the  flat  wall  case. 
The  composite  skin  friction  coefficient  is  as  much  as  56%  higher 
and  is  an  average  of  42%  higher  than  the  flat  wall  case.  The 
average  increases  in  Cf-  in  the  nonturbulent  and  turbulent  zone  are 
45%  and  23%,  respectively.  These  results  indicate  that  wall  cur¬ 
vature,  even  in  the  presence  of  strong  acceleration  and  high  FSTI, 
can  lead  to  a  significant  increase  in  momentum  transport  in  the 
boundary  layer.  The  Cf  values  from  the  mild  curvature  case  fall 
between  the  flat  wall  and  strong  curvature  results  as  expected.  The 
differences  between  the  flat  and  curved  wall  cases  shown  in  Figs. 
5-9  are  apparent  even  at  the  most  upstream  stations.  This  is  not 
unreasonable  since  G  =  5,  and  (rl  ds)0S=25  by  Station  1  of  the 
strong  curvature  case. 

Fluctuating  Velocity.  Figure  10  shows  u'  profiles  in  wall 
coordinates.  The  strong  curvature  composite  flow  results,  shown 
in  Fig.  10(a),  are  typical  of  transitional  and  turbulent  boundary 
layers.  The  peak  in  u' luT  occurs  at  y  +  =  15,  and  the  magnitude  of 
the  peak  is  between  2  and  2.3.  The  nonturbulent  and  turbulent 
zone  data  are  shown  in  Figs.  10(b)  and  10(c),  respectively.  In¬ 
cluded  in  these  figures  are  bars  showing  the  range  of  flat  wall  data 
from  Volino  et  al.  [19].  In  the  nonturbulent  zone,  the  concave  wall 
data  show  good  agreement  with  the  flat  wall  results.  In  the  turbu¬ 
lent  zone  the  data  from  all  stations  collapse,  particularly  near  the 
wall,  showing  self-similarity  throughout  the  transition  region.  In 
the  near  wall  region  u' tur  on  the  concave  wall  is  slightly  but 
consistently  reduced  compared  to  the  flat  wall  case.  This  indicates 
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Fig.  12  Fluctuating  wall-normal  velocity  profiles  in  wall  coor¬ 
dinates:  (a)  composite;  (b)  nonturbulent;  (c)  turbulent 


higher  skin  friction  relative  to  near  wall  turbulence  on  the  concave 
wall,  and  may  be  due  to  higher  mixing  in  the  outer  part  of  the 
boundary  layer.  In  the  outer  region  the  two  data  sets  show  good 
agreement  in  wall  coordinates.  Figure  11  shows  the  u'  profiles  at 
Station  5,  midway  through  transition  on  the  concave  wall.  The 
peak  in  u'  is  slightly  closer  to  the  wall  in  the  turbulent  zone,  but 
the  magnitudes  of  the  turbulent  and  nonturbulent  peaks  are  very 
nearly  equal.  The  difference  between  the  turbulent  and  nonturbu¬ 
lent  zones  is  much  reduced  from  the  flat  wall  case,  indicating  that 
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Fig.  13  Fluctuating  wall-normal  velocity  profile  for  Station  5 


something  in  the  concave  wall  case  is  promoting  fluctuations  in 
the  streamwise  direction  independent  of  turbulence. 

Figure  12  shows  the  v'  profiles.  Figure  12(a)  presents  the  com¬ 
posite  flow  results  for  the  strong  curvature  case.  The  profiles  are 
typical  for  a  high  FSTI  boundary  layer.  They  show  a  peak  v'  in 
the  near  wall  region  with  a  minimum  farther  from  the  wall  fol¬ 
lowed  by  a  rise  to  the  freestream  value.  The  minimum  in  U'  re¬ 
sults  as  the  freestream  turbulence  is  damped  by  the  wall.  The 
magnitude  of  v'/uT  rises  in  the  near  wall  region  as  transition 
proceeds.  As  was  also  observed  in  the  flat  wall  case  of  Volino 
et  al.  [19],  v1  is  lower  in  magnitude  and  shows  more  change 
through  the  transition  region  than  does  u' .  As  was  the  case  for 
u'luT,  the  turbulent  zone  results  (Fig.  12(c))  for  v’ tur  collapse 
for  all  of  the  stations.  The  major  difference  between  the  present 
results  and  the  flat  plate  results  is  the  high  value  of  v' lur  in  the 
turbulent  zone  in  the  outer  part  of  the  boundary  layer  and  the 
freestream  on  the  concave  wall.  Another  difference  is  the  reduc¬ 
tion  in  near  wall  v' luT  in  both  zones  relative  to  the  flat  wall  case. 
As  was  observed  for  the  turbulent  zone  u'/uT  profiles,  higher 
mixing  in  the  outer  part  of  the  concave  wall  boundary  layer  results 
in  higher  mean  velocity  gradients  at  the  wall,  raising  the  wall 
shear  stress  and  uT .  Higher  uT  reduces  both  v' /uT  and  u' luT  near 
the  wall. 

Figure  13  shows  the  U '  profiles  at  Station  5.  The  magnitude  of 
v '  in  the  turbulent  zone  is  almost  double  that  in  the  nonturbulent 
zone  across  the  entire  profile.  This  is  in  contrast  to  the  17'  profiles 
of  Fig.  1 1  that  show  the  two  zones  to  have  very  similar  magni¬ 
tudes.  The  same  difference  is  also  seen  in  the  flat  wall  case.  Volino 
[28]  found  that  much  of  the  unsteadiness  in  u  is  low-frequency 
unsteadiness  induced  by  the  freestream,  and  is  a  feature  of  both 
the  turbulent  and  nonturbulent  zones.  The  v'  fluctuations,  in  con¬ 
trast,  are  more  closely  related  to  near  wall  produced  turbulence 
and  eddy  transport  in  the  boundary  layer.  The  differences  between 
the  turbulent  and  nonturbulent  zones  explain  the  rise  in  near  wall 
v'  in  Fig.  12(a)  as  transition  proceeds.  Figure  13  also  shows 
differences  in  v '  between  the  concave  and  flat  wall  results.  The 
concave  wall  v'  profiles  for  both  zones  rise  more  rapidly  near  the 
wall,  reach  a  minimum  closer  to  the  wall,  and  are  higher  at  the 
edge  of  the  boundary  layer. 

Turbulent  Shear  Stress.  Turbulent  shear  stress  profiles  are 
shown  in  Fig.  14  in  wall  coordinates.  In  the  composite  profiles 
(Fig.  14(a)),  the  normalized  turbulent  shear  stress  increases 
through  transition.  In  the  nonturbulent  zone  the  turbulent  shear 
stress  is  low  and  comparable  to  the  flat  wall  case  results.  In  the 
turbulent  zone  (Fig.  14(c)),  there  is  a  significant  increase  in  the 
normalized  shear  stress  in  the  outer  part  of  the  boundary  layer 
compared  to  the  flat  wall  case.  These  outer  region  values  increase 
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Fig.  14  Reynolds  shear  stress  profiles  in  wall  coordinates:  (a) 
composite;  (to)  nonturbulent;  (c)  turbulent 


as  the  flow  moves  downstream,  presumably  due  to  the  increasing 
strength  of  curvature  as  the  boundary  layer  thickens.  Figure  15 
presents  the  turbulent  and  nonturbulent  —u'v'  profiles  for  Station 
5.  Similar  to  v'  and  in  contrast  to  u' ,  the  magnitude  of  —u'v1  is 
much  higher  in  the  turbulent  zone  than  the  nonturbulent.  Although 
—  u'v'  is  much  smaller  in  the  nonturbulent  zone,  it  is  nonzero. 
This  indicates  that  some  eddy  transport  of  momentum  occurs  even 
when  the  boundary  layer  is  nonturbulent.  The  figure  also  shows 
that  turbulent  zone  —u'v'  is  about  20%  higher  on  the  concave 


Fig.  15  Reynolds  shear  stress  profile  for  Station  5 


wall  than  on  the  flat  wall.  The  composite  flow  —u'v'  profile  is 
higher  on  the  concave  wall  due  to  it  being  further  through  transi¬ 
tion  than  the  flat  wall  profile  at  the  same  streamwise  location. 

Profiles  of  the  eddy  viscosity  are  presented  in  Fig.  16.  The 
composite  profiles  in  Fig.  16(a)  show  that  sM  increases  in  the 
streamwise  direction  as  the  transition  proceeds.  Comparison  of  the 
nonturbulent  and  turbulent  zone  profiles  in  Figs.  16 (to)  and  16(c) 
shows  that  the  eddy  viscosity  is  much  higher  in  the  turbulent  zone. 
While  the  difference  between  the  two  zones  is  clear,  Fig.  16(b) 
shows  that  there  is  still  significant  eddy  transport  in  the  nonturbu¬ 
lent  zone.  It  is  also  of  note  that  the  eddy  diffusivity  is  significantly 
increased  on  the  concave  wall  compared  to  the  flat  wall  in  both 
the  turbulent  and  nonturbulent  zones.  As  was  noted  for  the  turbu¬ 
lent  shear  stress  of  Fig.  14(c),  the  turbulent  zone  sM  of  Fig.  16(c) 
increases  at  the  downstream  stations,  due  to  the  rising  strength  of 
curvature  as  the  boundary  layer  grows. 

Streamwise  Vortices.  A  possible  explanation  for  the  differ¬ 
ences  between  the  concave  and  flat  wall  cases  might  be  the  pres¬ 
ence  of  Gortler  vortices  in  the  concave  wall  cases.  The  Gortler 
number  at  Station  3  of  the  strong  curvature  case  is  6,  indicating 
that  it  is  possible  that  vortices  could  form.  No  evidence  of  station¬ 
ary  Gortler  vortices  was  observed,  however.  Stationary  vortices 
would  result  in  spanwise  variation  in  velocity  profiles,  but  veloc¬ 
ity  profiles  measured  at  different  spanwise  location  and  the  same 
streamwise  location  were  indistinguishable.  Attempts  at  flow  vi¬ 
sualization  using  smoke  did  not  reveal  Gortler  vortices,  although 
the  high  FSTI  tended  to  scatter  the  smoke  quickly,  making  the 
flow  visualization  results  inconclusive.  If  vortices  were  present  it 
is  unlikely  that  they  would  remain  stationary.  In  a  low  FSTI  case 
vortex  location  would  be  fixed  by  some  small  upstream  distur¬ 
bance  such  as  an  upstream  screen  in  a  wind  tunnel.  In  the  present 
case,  the  turbulence  generating  grid  produces  large,  unsteady  fluc¬ 
tuations  which  might  be  expected  to  induce  nonstationary  vorti¬ 
ces.  Nonstationary  vortices  might  result  in  high  u'  away  from  the 
wall  as  the  upwash  and  downwash  locations  of  the  vortices  moved 
in  the  spanwise  direction.  Spectral  analysis  of  if'  in  the  present 
case  showed  a  broadband  peak  centered  at  30  Hz,  which  was  not 
present  in  the  flat  wall  case.  This  along  with  the  high  v'  in  the 
outer  region  (Fig.  12)  suggests  that  nonstationary  vortices  may  be 
present. 

Conclusions 

The  effect  of  concave  curvature  on  transitional  boundary  layers 
subject  to  high  freestream  turbulence  and  strong  favorable  pres¬ 
sure  gradients  has  been  documented.  Despite  the  strong  effects  of 
the  freestream  turbulence  and  acceleration  on  transition,  curvature 
still  has  a  significant  effect.  Curvature  causes  higher  momentum 
transport  in  the  outer  part  of  the  boundary  layer,  resulting  in  a 
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Fig.  16  Eddy  viscosity  profiles:  (a)  composite;  (b)  nonturbu- 
lent;  (c)  turbulent 


more  rapid  transition  to  turbulence  and  higher  skin  friction.  Con¬ 
ditional  sampling  shows  that  the  curvature  effect  is  present  in  both 
the  turbulent  and  non-turbulent  zones  of  the  transitional  flow.  In 
the  nonturbulent  zone  concave  curvature  produces  mean  and  fluc¬ 
tuating  velocity  profiles  that  appear  less  laminar-like  than  in  an 
equivalent  flat  wall  case.  In  the  turbulent  zone  curvature  produces 
higher  fluctuating  velocity  and  turbulent  shear  stress  in  the  outer 
part  of  the  boundary  layer.  Nonstationary  Gortler  vortices  are  of¬ 
fered  as  a  possible  mechanism  for  the  curvature  effect  on  the 
boundary  layer.  Although  some  results  which  suggests  the  pres¬ 


ence  of  vortices  are  presented,  no  direct  evidence  of  their  exis¬ 
tence  was  found.  Further  study  is  needed  to  confirm  the  mecha¬ 
nism  for  the  concave  curvature  effect. 
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Nomenclature 


Cf  = 

FSTI  = 

f(Vpk)  = 
G  = 
H  = 
K  = 
n  = 
Re*  = 
Re,  = 
r  = 
U  = 
U_w  = 
t/„  = 
u  = 

u+  = 


uT  = 

—  u'v'  — 

—  u'v'  — 

v'  = 

X  = 

y  = 

y+ : 
^99.5  — 

S*  = 
em  = 

y  = 

y Pk  = 

V  = 

p  = 
0  = 
cr  = 


Subscripts 


tw / (pU\j2) ,  skin  friction  coefficient 
freestream  turbulence  intensity 
function  of  intermittency,  Eq.  (1) 

Re e(8/r)05,  Gortler  number 
S*/ 0,  shape  factor 

(v/U2,)(dUw/dx),  acceleration  parameter 
dimensionless  turbulent  spot  production  rate 
Uwxlv,  Reynolds  number 
momentum  thickness  Reynolds  number 
radius  of  curvature  of  test  wall 
time-averaged  local  streamwise  velocity 
local  freestream  velocity  extrapolated  to  the  wall 
average  freestream  velocity  in  transition  region 
instantaneous  streamwise  velocity 
U/uT,  local  mean  streamwise  velocity  in  wall  coordi¬ 
nates 

rms  streamwise  fluctuating  velocity,  yju'2 
\J  tw  /  p ,  friction  velocity 
instantaneous  turbulent  shear  stress 
time-averaged  turbulent  shear  stress 

rms  cross-stream  fluctuating  velocity,  \jv  '2 
streamwise  coordinate,  distance  from  leading  edge 
cross-stream  coordinate,  distance  from  wall 
yuT/v,  distance  from  wall-in-wall  coordinates 
99.5%  boundary  layer  thickness 
displacement  thickness 
—  u'v' /(dU/dy),  eddy  viscosity 
time-averaged  intermittency,  fraction  of  time  flow  is 
turbulent 

peak  intermittency  in  profile 
kinematic  viscosity 
density 

momentum  thickness 

turbulent  spot  propagation  parameter 

wall  shear  stress 


■j  =  transition  start 
e  =  transition  end 
co  =  local  freestream  condition 
w  =  extrapolated  to  the  wall 
NT  =  nonturbulent  zone 
TURB  =  turbulent  zone 
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